Finite size effects of helical edge states in HgTe/CdTe quantum wells 
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The solutions for the helical edge states for an effective continuum model for the quantum spin 
Hall effect in HgTe/CdTe quantum wells are presented. For a sample of a large size, the solution 
gives the linear dispersion for the edge states. However, in a finite strip geometry, the edge states 
at two sides will couple with each other, which leads to a finite energy gap in the spectra. The gap 
decays in an exponential law of the width of sample. The magnetic field dependence of the edge 
states illustrates the difference of the edge states from those of a conventional quantum Hall strip 
of two-dimensional electron gas. 

PACS numbers: 73.43.-f, 72.25.Dc, 85.75.-d 



Recent discovery of quantum spin Hall (QSH) effect 
brings the Hall family a new member, an insulator with 
topological properties of electron bands distinct from the 
conventional ones Kane and Mele 0] proposed that 
newly synthesized graphene may exhibit a QSH effect. 
This novel state consists of the helical edge states whose 
dispersions locate inside the bulk insulating gap. They 
indicated that the Z2 topological number distinguishes 
between the QSH and usual insulating phases. Benevig 
et al [c| realized that the HgTe/CdTe quantum well is 
a potential candidate for observing this effect based on 
the fact that the electron bands changes from the normal 
to an " inverted" type at a critical thickness of the quan- 
tum well, which was confirmed experimentally within one 
year.Q Several other candidates were also proposed, such 
as GaAs with shear strain [B| and a multilayer Bi thin 
filmjfl. Up to now, considerable efforts have been done 
to explore the properties of the QSH effect and topolog- 
ical insulators 7]. One of the striking properties is the 
crossing linear dispersions of the two edge states in which 
the electric currents with different spins flow in opposite 
directions. Very recently, it was reported that the mas- 
sive Dirac particles exist in the bulk of Bio.gSbo.i, which 
is a hall mark of higher dimensional quantum spin Hall 
insulator. Q 

The theory of QSH effect in HgTe/CdTe was based on 
an effective 4-band model of the HgTe/CdTe quantum 
well that participates in the inversion crossing of elec- 
tron and hole bands derived by Benervig, Hughes and 
Zhang [3] . Until present almost all works are based on nu- 
merical solutions of the model in a tight-binding method 
It was argued that the QSH insulator state has 
gapless and linear dispersions locating inside the bulk in- 
sulating gap. In the present Letter, we present the solu- 
tion for the edge states of the 4-band model. In a finite 
strip geometry, linear dispersions for four edge states are 
reproduced for a sufficient large sample comparing with 



the space distribution scale of the edge states. However, 
when the size of sample is comparable with the size scale 
of the edge state, the edge states near the crossing point 
will couple with each other and open a finite energy gap. 
The gap decays in an exponential law with the size of the 
sample. This demonstrates that the edge states of QSH 
effect are quite different from the edge states of quantum 
Hall (QH) effect. 

Here we start from the effective 4-band model for 
HgTe/CdTe quantum wells [3] 
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where H (fc) = £fel 2 + d a (k) a a , with I 2 being a 2 x 2 unit 
matrix, <r a the Pauli matrices, e& = C—D (fc 2 , + ky) , d 1 — 
Ak x , d 2 = Ak y , and d 3 = M (k) = M — B (fc 2 



A, B, C, D, and M are material specific parameters 
that are functions of the thickness of the quantum well. 
The model is derived from the Kane model near the T 
point in a heterojunction of HgTe/CdTe, and its valid- 
ity is limited for small values of k XiV . The most strik- 
ing property of this system is that the mass or gap pa- 
rameter M changes sign when the thickness d of the 
quantum well is varied through a critical thickness d c 
(= 6.3 nm) associating with the transition of electronic 
band structure from a normal to an " inverted" type 11 1 . 
We will now solve this continuum model {T]) in a finite 
strip geometry of the width L with the periodic bound- 
ary condition in the x direction and an open bound- 
ary condition in the y direction. In this case, k x is 
a good quantum number but k y is replaced by using 
the Peierls substitution k y — —id/dy = —id y . The 
Hamiltonian ([I]) is block-diagonal, and the upper iff 
(= H (fc)) and lower H{ (= H* (—k)) blocks describe 
the states of spin-up and spin-down, respectively. The 
eigenvalue problem of the upper and lower blocks can 
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be solved separately, i.e., H^^ = and H^i = 

with *i(k x ,y) = e ik " s {if, 1 (k x ,y),il> 3 (k x ,y)) T 
and *| (k x ,y) = e ik * x (i/) 3 {k x , y) , ip 4 (k x ,y)f where T 
means "transpose". Because the lower block of the 
Hamiltonian is the time reversal of the upper block of 
the Hamiltonian, the solution (k x ,y) = OvI/j (k x ,y), 
where = —i<j y K is a "time-reversal" operator and K 
stands for complex conjugation. Thus we can only focus 
on the solution for the upper block of this Hamiltonian. 

The set of the eigenvalue equations for the upper block 
is expressed as 

[M-B + {kl-dl)]-4> 1 +A{k w -d y )^ 2 = £Vi,(2) 
A (k x + d v ) ^ - [M - B_ (kl - d 2 y )} V- 2 = £V 2 ,(3) 

with B± — B±D. From eqs. ^ and ([3]), using the trial 
function ^12— eX% '■> t ne characteristic equation gives four 
roots ±Ai and ±A2, 



where F 



K,2 = k 2 x +F± ^F 2 -{M 2 -E 2 )/B+B-. (4) 

A 2 -2(MB+ED) 



With the boundary condi- 



2B + B_ 

tions of (k x ,y = ±L/2) = 0, we have an analytical 
expression for the wave function ^| 



ipi = (k x ,y) + c-f- 
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d+f+ (k x ,y) + d-f- 



(k x ,y) , 
(kx,y) , 



with 



f+(k x ,y) = 
f- {k x ,y) = 



cosh(Aij/) cosh(A2y) 
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sinh (Aij/) 
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sinh (X 2 y) 
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(5) 
(6) 
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The non-trivial solution for the coefficients c± and d± in 
the wave functions leads to a secular equation 
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where a(E) = E - M 



B+k 2 x 



E-M + B+k 2 x - B+Xj. Eqs. g| and 



B+X{, and P{E) = 
([9]) determine the 
energy spectra of the upper block Hp 

Let us first consider the general properties of the solu- 
tion for Ai i2 . In a large L limit, a purely imaginary A = 
iky is always the solutions of the equation, and gives two 



branches of spectra, E± = e k ± J (M - Bk 2 ) +A 2 k 2 . 
These are the bulk spectra and revised by the A-term 
slightly, and corresponding solutions span in the whole 
space. Except for these imaginary solutions, there also 
exist real solutions when 



A 2 /B + B_ > max{2M/B,4M/B}. 



(10) 



In the large L limit, tanh (Ai^i/2) = 1, and Eq. (j9| 
gives 

E± = M - B + XiX 2 ± B+(Ai + X 2 )k x - B+k 2 x . (11) 
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FIG. 1: (Color online) (a) Energy spectra of edge states for 
L — 200 nm (blue solid lines) and L = 1000 nm (red dashed 
lines) by solving eqs. Q and ([9]) for the HgTe/CdTe quan- 
tum well thickness d = 7.0 nm. Fhe parameters are given in 
the text. As a comparison, numerical results of tight binding 
approximation, which was used in Ref.|9| for the same param- 
eters are also plotted as black squares for L = 200 nm and 
black dots for L — 1000 nm. (b) Fhe width dependence of 
the energy gap A for a finite width strip of HgTe/CdTe quan- 
tum well of thickness d = 7.0 nm. Specifically, A(L = 200 
nm) = 0.4509 meV, A(500 nm) = 1.6 x 10~ 3 meV, and 
A(L = 1000 nm) = 1.41 x 10" 7 meV. 



As the Ai_2 approaches to constant near k x = 0, 
E± (k x ) ~ -MD/B ± Ay / B + B^/B 2 k x . For real roots 
A1.2, the function f±(y) are distributed dominantly near 
the edge (y = ±L/2) in the scale of X^\. Fhis result 
is consistent with those by means of the tight binding 
approximation [9(. 

In the following we shall concentrate on the real solu- 
tion of A with a finite L, i.e. the solutions for the edge 
states. For real A and finite L, the right hand side of Eq. 
© is always greater than 2. If Xi 2 L >> 1, it is approxi- 
mately 2 + 4e~ 2A2i (assuming Ai >> A 2 ). From Eq. ([9]), 
it is found that a finite energy gap A = E+ — E- opens 
at k x = 



±\AB+B-M\ 



-X 2 L 



y/B 3 (A 2 B - AB+B-M) 



(12) 



which decays in an exponential law of L. This is the main 
consequence in the present work. 

In general cases, we have numerical solution of the 
equations. As a concrete example, we adopt the param- 
eters for the inverted HgTe/CdTe quantum well of thick- 
ness d = 7.0 nm from the reference [10] for all numerical 
calculations in the present Letter: A = 364.5 meV nm, 
B = -686 meV nm 2 , M = -10 meV, D = -512 meV 
nm 2 . For L = 1000 nm, one has Ai = 0.7797 nm' 1 and 
A2 = 0.0187 nm -1 at k x = 0. The energy gap is very 
tiny, A = E + — £L = 1.41 x 10~ 7 meV. However, for 
L = 200 nm at k x = 0, the gap A = 0.4509 meV, which 
becomes large enough to be measurable in experiments. 
We plot the energy spectra for the edge states of several 
sizes in Fig. 1(a). However, for a narrow width L (e.g., 
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FIG. 2: (Color online) The density distribution of the two 
edge states (k x ,y) for L = 200 nm. (a) The red solid 
line corresponds to |*1/|+ (k x , y)\ 2 and the blue dotted line 
to * T+ (-k x ,y)\ 2 at k x = 0.01 nm" 1 ; (b) The red solid 
line corresponds to |\P T _ (k x ,y)\ 2 and the blue dotted line to 
|* T _ (~k x ,y)\ 2 at k x = -0.01 nm -1 ; (c) and (d) for k x = 



L = 200 nm), there is the parabolic-like spectrum near 
k x = 0. The size-dependence of the energy gap is plotted 
in Fig. 1(b). 

The corresponding wave functions of eigenvalues 
E±{k x ) yields 



where 
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c+e ik " (f ++1 ij-,r,tf-+livtf + ) C13) 



B + (A 2 - A 2 ) I A 



A 2 L 



Ai coth ^ - A 2 coth „ 
B + (A 2 - Xt) I A 



Ai tanh ^ - A 2 tanh 
B + - x i) k x r]i/V2 



E=E± 



E=E± 



/3Ai tanh ^ - aA 2 tanh ^ 
5+ (A 2 - Aa) KVzhx 



/?Ai coth - aA 2 coth ^ 



E=E + 



E=E- 



and c± are normalization constants. The solutions can 
be simplified in the limit of large L as r/f = 4 = Tj = 
B+(Ai + A 2 )/A and 7^ = -7^ = -sgn(k x ). The other 
two solutions for the lower block can be produced by 
means of the time reversal operation, ^ i± — Q^]±, and 
the spectra are degenerate with those of the upper block. 

According to the present analytic solutions of wave 
functions, the density distribution of the functions are 
mainly determined by the two length scales, Xi\. In the 
example in Fig. 1, we notice that A^ 1 >> X^ . For a 



large size of the sample, the density of the wave func- 
tion increases in an exponential law in the scale of A^ 1 
and then decays exponentially in a scale of X^ 1 near the 
boundaries, which is consistent with the work by Konig et 
al.\ld^ The wave function almost vanishes far away from 
the boundaries if the width of the sample is much larger 
than X% . As an example, the density distributions of 
^T± (fee, y) for L = 200 nm are plotted for demonstration 
in Fig. 2 where A^ 1 — 55.9 nm and 51.8 nm at k x = 0. 
The states of <Fj + (k x ,y) and *&-[■_ (— k x , y) (k x > 0) have 
the same spin (cx (1, —rj) T in the large L limit) and the 
positive velocity, v x (= +A\/B + B-/B 2 ) > when k x 
is far away from k x = and the density distribution 
is located at one side while the states of ^-\+ (— k x , y) 
and 'Fj- (+k x ,y) (k x > 0) have another spin (oc (l 7 r]) T 
in the large L limit) and a negative velocity, — v x < 0, 
and are distributed on the other side. From the solution 
we found that \I/|± (k x ,y) and ^]-± (—k x ,y) couple near 
k x = due to the finite size effect. Consequently, the 
states with different spins mix together and the densities 
of the wave functions 'Fji (k x = 0,y) are symmetrically 
distributed at the two sides. This fact is consistent with 
the opening of an energy gap in the spectra at k x = 0. 
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FIG. 3: (Color online) The variations of conductance G via 
the chemical potential fi inside the bulk insulting gap for a 
HgTe/CdTe quantum well of L — 200 nm at temperatures 
T = 30 mK and 1.8 K, respectively. Note that the energy 
zero point is shifted to the center of the gap A (= 0.4509 
meV). 



The charge conductance of a QSH phase in a strip was 
predicted theoretically to be 2e 2 /h due to the presence 
of two spin-resolved conducting channels at the edges of 
the strip, which was observed experimentally in two sam- 
ples with sizes of (1.0 x 1.0) ^m 2 and (1.0 x 0.5) ^m 2 , 
respectively. [1] The finite size effect will modify the con- 
ductance of the QSH phase. Following the Landauer- 
Biittiker formula, 
form, 



12j the charge conductance has the 



2e 2 

0(A)-* 
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G (A) — > 2e 2 /h at low temperatures only when the chem- 
ical potential locates out of the gap A(L). Below the tem- 
perature of k B T* = A, a dip will be obviously exhibited. 
We plot the temperature dependence of the conductance 
in Fig. 3. In the experiment by Konig et aZ[4], the small- 
est sample has L = 500 nm, and the measurement was 
performed at 30 mK. The calculated A = 1.6 x 10~ 3 meV 
(19 mK), which is already comparable with the experi- 
ment temperature. For a smaller example, the energy gap 
of L = 200 nm is about A = 0.4509 meV, and the corre- 
sponding temperature is enhanced to T* — 5.22 K. Ex- 
perimentally, the data for two smaller samples are close 
to the value of 2e 2 /h while the data for the wider sam- 
ples obviously deviate from the value. From the present 
exact solution, it is believed that the value will be also 
modified at lower temperatures for even smaller samples. 
By the way, it is also worth noting that the gap is also 
highly sensitive to the thickness of quantum well, as all 
parameters in Eq. (fTJ) are functions of the thickness. 
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FIG. 4: (Color online) The energy dispersion of the edge 
states in a weak field of B — 0.04T. The blue lines are for 
the upper block and red lines are for the lower block. The 
back lines is for the bulk spectra. The width of sample is 
L = 200nm. 

Finally, our solution shows that the QSH edge states 
is quite different from the edge states of a conventional 
QH strip. For a QH strip with translational symmetry 
along the strip, the states are classified with the momen- 
tum k x . The edge states at the two sides have different 
fe-lIU So they do not mix together even when the two 
states overlap in the space if there is no other scatter- 
ers or interactions in between two edges. In the case of 
a QSH strip, the edge states at the two sides have the 
same k x . Near the anti-crossing points, the two states 
has the nearly equal energy and momentum. So they 
can couple together to generate an energy gap when their 
energies becomes closer and the wave functions have over- 
laps in a finite space. The magentic field dependence of 
the QSH edge states also reflects this peculiar property. 
Consider the sample is subjected to a weak perpendic- 
ular magnetic field B z . Using the Peierls substitution, 
k x —> k x — eA x /h in Eq.(l) by taking the gauge, A x = 



—B z y (for \y\ < L/2) and A y = in order to keep k x a 
good quantum number. As the wave functions of the edge 
states decay exponentially, the expectation value of A x in 
the two edge states is proportional to L approximately 
for L » A7.2- The energy will be shifted by AE = 
+gpL B B z (g « m e v x [L - Ar 1 - A7 2 - 2(Ai + A 2 ) _1 ] /h) 
near k x = 0. Thus the energy spectra of the two edge 
states of the upper block will shift downward or upward 
E{k x ) « ±v x hk x +g^, B B z for a large L. By increasing the 
magnetic field B z the anti-crossing point of energy spec- 
tra is eventually moved out of the bulk insulating gap, 
and the spectra will not crossing in momentum between 
the gap. However, another two branches of the spectra 
of the lower block will move in an opposite direction, and 
the two sets of the spectra may cross in momentum inside 
the insulating gap. A more detailed calculation is to use 
the exact solutions of Eqs. (13) and (14) at B z = as a 
basis to truncate the model of Eq.(l) in a B z field to an 
effective one. Numerical results are plotted in Fig. 4. The 
energy shift is very sensitive to the magnetic field and the 
width of the sample. As for the conductance, the value 
of 2e 2 /h will recover near the crossing points. Thus the 
magnetoresistance is very sensitive to a tiny field, which 
might have potential application for a sensitive detection. 
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